Partitioning Regular Polygons into Circular Pieces II: 

Nonconvex Partitions 



Mirela Damian* Joseph O'Rourke^ 



Abstract 

We explore optimal circular nonconvex partitions of regular fc-gons. The circularity of a polygon 
is measured by its aspect ratio: the ratio of the diameters of the smallest circumscribing circle to the 
largest inscribed disk. An optimal circular partition minimizes the maximum ratio over all pieces in 
the partition. We show that the equilateral triangle has an optimal 4-piece nonconvex partition, the 
square an optimal 13-piece nonconvex partition, and the pentagon has an optimal nonconvex partition 
with more than 20 thousand pieces. For hexagons and beyond, we provide a general algorithm that 
approaches optimality, but does not achieve it. 

1 Introduction 

In |DO03| we explored partitioning regular fc-gons into "circular" convex pieces. Circularity of a polygon is 
measured by the aspect ratio: the ratio of the diameters of the smallest circumscribing circle to the largest 
inscribed disk. We seek partitions with aspect ratio close to 1, ideally the optimal ratio. Although we start 
with regular polygons, most of the machinery developed extends to arbitrary polygons. 

For convex pieces, we showed in jDU03j that optimality can be achieved for an equilateral triangle only 
by an infinite partition, and that for all fc > 5, the 1-piece partition is optimal. We left the difficult case of 
a square unsettled, narrowing the optimal ratio to a small range. Here we turn our attention to partitions 
that permit the pieces to be nonconvex. The results are cleanest if we do not demand that the pieces 
be polygonal, but rather permit curved sides to the pieces. The results change dramatically compared to 
the convex case. The equilateral triangle has an optimal 4-piece partition, the square an optimal 13-piece 
partition, the pcntagan an optimal partition with more than 20 thousand pieces. For hexagons and beyond, 
we provide a general algorithm that approaches optimality, but does not achieve it. 

1.1 Notation 

A nonconvex partition of a polygon P is a collection of sets Si satisfying 

1. Each Si C P. 

2. UA = p. 

3. The sets have pairwise disjoint interiors. 

These conditions alone are too broad for our purposes, as there are no constraints placed on the pieces. It is 
reasonable to demand that each set be connected, but even this is too broad. The most natural constraint 
for our purposes is to require the interior of each piece to be connected: 
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4. The interior of each Si is connected. 

The aspect ratio of a piece is the ratio of the radius of the smallest circumcircle to the radius of the largest 
inscribed disk. Aspect ratios will be denoted by symbol 7, modified by subscripts and superscripts as 
appropriate: 71 (-P) is the one-piece 7; "f{P) is the maximum of all the 71 (S'^) for all pieces Si in a partition 
of P; J*{P) is the minimum 7(P) over all nonconvex partitions of P. Our goal is to find "y*{P) for the 
regular fc-gons P. Both the partition and the argument "(P)" will often be dropped when clear from the 
context. 

Throughout we consider all disks to be closed sets, including the points on their bounding circle. Disks 
will be denoted either by symbols i — 1,2,3, ... ,n; the subscript will indicate the disk bound by an 
inscribed/in-disk, and 1 will indicate the circumscribed/out-circle. 

1.2 Table of Results 

Our results are summarized in Table [L^ 



Polygon 


71 


79 


nonconvex, nonpolygonal 


7* 


fc* 


Equilateral Triangle 


2.00000 


1.50000 


7e 


4 


Square 


1.41421 


1.20711 


79 


13 


Regular Pentagon 


1.23607 


1.11803 


79 


< 20476 


Regular Hexagon 


1.1547 


1.07735 


1.10418 


finite 


Regular Heptagon 


1.10992 


1.05496 


1.08382 


finite 


Regular Octagon 


1.08239 


1.0412 


-yi = 1.08239 


finite 


Regular fc-gon 


1/ cos(7r/fc) 


l + csc(e/2) 
2 


< 7f = 1.08239 


finite 



Table 1: Table of Results on Regular Polygons. 71: one-piece partition; 7]"': one piece ratio 71 for regular 
fc-gon; 7g: single-angle lower bound; 9: angle at corner; 7*: optimal partition; k*: number of pieces in 
optimal partition. 



Here we use to denote the "one-angle lower bound" , a lower bound derived from one angle 9 of the 
polygon, ignoring all else. This presents a trivial lower bound on the aspect ratio of any partition. 

2 Preliminary lemmas 

We recall two simple lemmas used in Table [T!^ proved in |UO03| . 
Lemma 1 (Regular Polygon). The aspect ratio 71 of a regular k-gon is 

1 

cos{n/k) 



Lemma 2 (One- Angle Lower Bound). // a polygon P contains a convex vertex of internal angle 
then the aspect ratio of a partition of P is no smaller than jg, with 

1 -hcsc(6i/2) 
7e = 7: 
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3 Equilateral Triangle 

An equilateral triangle has 71 = 2. The lower bound provided by Lemma |21 is 79 = 1.5 (see Table [TT^ . 
Figure^shows a partition with 4 pieces that achieves 79, and is therefore optimal. This partition has three 
convex corner pieces and one nonconvex central piece. 




Figure 1: Optimal partition of an equilateral triangle (4 pieces). Inscribed and circumscribed circles are 
shown. 

4 Square 

A square has 71 = ^2 « 1.41421. The lower bound provided by LemmaElis 79 = (1 + v^)/2 « 1.20711 
(see Table IT!^ . Figure shows a partition with 13 pieces that achieves 76/, and is therefore optimal. The 




(a) (b) 
Figure 2: (a) Optimal partition of a square (13 pieces) (b) Magnified view of one square corner. 

partition contains one large central nonconvex piece, four convex corner pieces and one nonconvex piece to 
each side of each corner piece, for a total of 13 pieces. 

As k increases, decreases and it becomes increasingly difficult to partition a fc-gon into pieces with 
optimal ratio. As hinted in the square partition, it becomes essential to be able to cover small gaps along 
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the interior of edges. Even for the pentagon, a less ad hoc procedure is needed. In the next section, we 
devise a general algorithm that covers a subsegment of an edge with pieces with ratio close to optimal. This 
will permit us to make progress for fc > 4. 

5 Covering an edge segment 

Let S be an edge segment tangent to two disks Dq and Dq at its endpoints. A covering of S is two collection 
of disks, Dq and Dl, i = 0,1,2,3, .. . with four properties: 

1. Each disk Dq is tangent to S 

2. The disks Dq have pairwise disjoint interiors: J{Dq) D J{Dq) ~ for i ^ j. 

3. The disks D\ collectively cover S: U^D\ D S. 

4. Each Dq is inside the corresponding D\: Dq C D\ for all i. 

For a given edge segment S, our goal is to find a covering of S of optimal ratio. In the following we present 
an algorithm that finds a covering of S of ratio close to the optimal. 

5.1 Algorithm (Edge Cover) 

The algorithm presented in this section takes as input: 

(a) An edge segment S — [aa,ai] 

(b) Disks Dq and Dq tangent to each other and to S at points oq and ai, respectively 

(c) Corresponding outcirclcs Di D Dq and Dl D Dq 

(d) The desired ratio factor 7 > 1 

and seeks to extend the sets {Dq, Dq} and {Di, Dl} to a covering of S of ratio 7, if one exists. 

For i £ {0, 1}, let ai be the point where Dq touches S and bi the point where D\ intersects S. We 
start by growing the largest possible indisk Dq that touches the uncovered segment piece at midpoint 
(i2 = {bo + bi)/2. Clearly, Dq can only grow until it touches either of the two adjacent indisks, Dq or Dq. 
We will show later that Dq hits the smaller of Dq and Dq first (see Figure Next we inflate Dq by 7 to 
obtain Df and displace D^ vertically downwards until its topmost point touches the topmost point of Dq, 
so as to capture as much of the uncovered edge segment as possible. 

If Dl covers the entire triangular gap (as in Figure (JJa), we are finished. Otherwise, recurse on the at 
most two new edge segments created: [ao, 02] and [02, oi]. Note that the uncovered gaps of these two edge 
segments are identical and therefore their coverings will be identical. 

5.2 Analysis 

Without loss of generality, we assume that Dq is at least as large as Dg. For analysis convenience, consider 
a coordinate system with the origin where intersects the horizontal edge, as in Figure 0] At a certain 
stage of the algorithm, all uncovered gaps in the original edge segment are symmetric and will be covered in 
the same way. In our analysis, we focus on the uncovered gap adjacent to the origin; henceforth, the term 
gap will refer to the leftmost uncovered gap of the edge segment, with leftmost understood. 

Let Dq be the indisk on the right side of the gap in iteration step n; Dq always remains to the left side 
of the gap. Refer to Figure^] For any n, let r„ denote the radius of Dq and let a„ be the point where Dq 
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(a) 



(b) 



Figure 3: Algorithm (a) Iterative step: Dq centered on the midpoint a2 of [bo-^i] (b) Termination: Df 
covers the gap. 

touches the x-axis. We define a useful quantity 5n to represent the distance from a„ to where D" intersects 
the X-axis: (5„ = r„ 1/72 — (2 — 7)^, or equivalently 



Sn = 2r„v/7^ (1) 

In iteration step n + 1, the algorithm grows the indisk D^^^ tangent to the uncovered gap [0,a„ — (5„] at 
its midpoint a„+i = (a„ — 5n)/2, until it hits either or Dq. Using (5„ from this is 




Proof: We determine r„+i from the tangency requirement (a„ — a„+i)^ + (r„ — r„+i)^ — (r„ + r„+i)^, or 
equivalently 

(a„-a„+i)2 

^^"+1 = Tv ' 

and show that i'g^"'^ and D\ are disjoint: 

(a„+i - ao)^ + (ro - r„+i)^ > (rg + r„+i)^ 
Substituting the expression for r„+i from © yields 



Note that ro > r„, since ro > r\ and r„ decreases as n increases. Also from Q we have ao = — 27'o-\/7 ^ 1- 
This together with (|21) renders the inequality above true. □ 
Our goal is to find the optimal 7 for which the algorithm terminates in a finite number of steps. 
This involves solving the coupled recurrence relations |2Il and Q and imposing the termination condition 
a„+i — (5„+i < 0, which ensures that the edge segment is completely covered in iteration step n + 1. 
Substituting 5„+i from (Q) yields 

which together with ||2Jl and Q leads to an system of recurrent relations with two variables. Next we show 
how to reduce these recurrence relations to only one recurrence relation in one variable, which is easily 
solvable. 

5.2.1 Rescaling the gap 

The leftmost segment gap we wish to cover is always bounded to the left by _Dq, whose position remains 
unchanged. This suggests a simple way to simplify the coupled recurrence relations Q and rescale Dg 
at the end of the iteration step n, so as to ensure r„ = 1 at the start of the iteration step n + 1. Initially, 
we scale the disk Dq and set — r\ j rx — \ and 

a, = - (4) 

Let a„ and r„ denote the scaled variables at the end of iteration step n, with r„ = 1. Based on (O and (j^J, 
we determine in iteration step n + 1 



^-VW (5) 

(6) 



' (a,! ^ '^n+\) 



Rescale D^'^^ to ensure = 1. Thus, < — r^^j^ / r„^i = 1 and 

a„+i < -, (7) 

Substituting in Q the expression for r„^2 from |(HJ yields one recurrence relation for a„ of the form 

aL+i = Fisl^) (8) 

with 



n-) = (9) 



2a; -4^7^ 

Lemma 0] establishes the relationship between the scaled a„ and its unsealed correspondent a„ 
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Lemma 4 For each n, = an / r„ at the end of iteration step n. 

Proof: The proof is by induction on n. The base case is n = 1, which is clearly true from (01 . Assume 
a„ = a„ / r„ for any n < s, for some s > 0. Now we show that a^j^-^ — Ogj^i / r^+i. We use the induction 
hypothesis = as/rs in (jHl to obtain ag_^i — F{as/rs). From ^ and ©, we get 



a^+i 4:rs{as/2 - r^^-f - 1) 



's+l 



[os - Os+lY 



Substituting again Os+i from ((21 in the expression above yields ag^i/rs+i — F{as/rs) = a^+i, which proves 



the lemma. 



□ 



5.2.2 Computing optimal 7 

The edge cover algorithm terminates when D" covers all points of the uncovered gap, i.e, a„ < (5„. From 
and the fact that r„ = 1, we derive the stopping condition 



(10) 



Our goal is to determine the optimal 7 for which inequality (|10() is satisfied for some finite n. Clearly, we 
want a„ to move down to 6 , getting closer to 5 with each iteration step; that is, a^j^-^ < a„ for all n. 
However, we show that this does not happen for any 7 and any edge segment [0, Ui]: 




Figure 5: Fixed point. 



Theorem 5 The algorithm terminates in a finite number of steps only if one of the following is true: 
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(a) 7 > 7* = 1.11340 

(b) 7 < 7* and F{ai) < and F (qi) > 1 

Proof: The proof consists of three parts. First we show that the equation F{x) — x has two positive roots 
a\> a\> 5 . Next we prove that the iteration procedure a„_,_]^ = F(a^ converges to a*, unless one of the 
two conditions (a) and (b) stated above is met. The imphcation of this is that the edge cover algorithm 
gets stuck at a\ and fails to make any further progress towards 6 ; hence, it never stops. Finally, we show 
that under either of the two conditions stated in the theorem, the algorithm terminates in a finite number 
of steps. 

Using @, we reduce x = F{x) to a cubic equation 



+ 4^7 - 1.^^ + 4(7 - 3)a: + 16^7 -1 = (11) 
which can be solved by use of Cardano's method jGul97| . Solving for x involves the determinant 

64 

A = --(47^- 797 + 83) (12) 
This quadratic polynomial has one root of interest 

7* = ^^-f^ = 1.11340 (13) 
8 

and a second root outside the domain of interest. We omit to show here the complicated expressions for 
the roots of equation lll|l . Figure [S] shows with solid curves how these roots vary with 7. The dashed line 
in Figure [5I shows the finishing point 8 . Note that for any 7 S [1,7*], the equation x = F(x) has three real 
roots: two positive roots a| > 02 > (5 , and one negative root. Figure |H| shows a magnified view of the two 
positive roots in the vicinity of 7* = 1.11340. 

We now show that for any 7 G [1,7*], the iteration procedure a„+i = i^(a„) converges to a|, unless 
condition (b) of the theorem holds. From the three regions delimited by the contours of the two positive 
roots al and in Figure observe the following: 

(a) if F{aj^) < Oi', then lies in region / above the curve F{al) — a|; therefore, F(a„) < a„ for some 

> and all n < N. Also note that F (x) < 1 in a neighborhood containing both al and a^, which 
guarantees that a„ converges to al < 

(b) if F{ai) > ai and F [a-^) < 1, then ai lies in region II delimited by the contours of the two positive 
roots; therefore, F{a„) > a„ for some > and all n < N. Again, since F (x) < 1 in a neighborhood 
containing both a* and a^, a„ converges to > a^. 

(c) if F{ai) < ai and F {a^) > 1, then lies in region /// below the curve -^(02) = flj. Hence, 
F(a„) < a„ for all n and therefore a„ reaches 5 in a countable number of steps. Also note that the 
same is true for any 7 > 7* (region IV in Figure E}. 

Finally, we show that if the algorithm terminates, then a„ reaches S in a finite number of steps. In other 
words, there exists a constant e > such that 

0,1+1 < a„ — £ 

is satisfied for any iteration step n. This is equivalent to 

F{x) <x-e (14) 

An analysis similar to the one of equation (|ll|l shows that there exists £ > that satisfies ()14|l for all x. 
This ensures that the number n of iteration steps is bounded above by {ai — 6 )/e. □ 
Figure shows the number of iteration steps it takes to cover a segment tangent on its endpoints to two 
unit radius disks tangent to each other. Note that for any 7 > 1.126, the edge segment can be covered in 
one step only; for any 1.116 < g < 1.126, the edge segment can be covered in two steps; and so on. As 7 
approaches the critical value 7*, the number of steps increases exponentially. 
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Figure 6: Fixed point. 
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Figure 7: Edge cover ratio 7 versus number of iterations. 
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5.3 Triangular gap partition 



Lemma 6 Let Dq and Dq be two disks tangent to each other and to an edge segment S at its endpoints. If 
Dq covers the intersection point between and D^, then a covering produced by the Edge Cover algorithm 
for S covers all points of the triangular gap delimited by S, Dq and Dq . 

Proof: Let t„ be the intersection point between ZJj and I?" closer to S in iteration step n (see Figure |H|- 
Note that i„ is the apex of the triangular gap left uncovered in iteration step n, which we attempt to cover 
in iteration step n+ 1. If for any n, Dq~^^ covers i„, then clearly the circles Dl collectively cover all points 
of the original triangular gap. 




Figure 8: Dq covers tn for all n > 1. 

As discussed earlier, the Edge Cover algoritm terminates only if a„^]^ < a„ for all n. Using Lemma 0] 
this is equivalent to 

^ rv^ (15) 

Q-n rn 

This tells us that a„ decreases at a faster rate than r„ with increasing n. Let qn be the intersection point 
between D" and the line L that passes through <i and origin, with qi = ti. Refer to Figure|Hl An implication 
of (|15|l is that that qn moves lower inside Dq with increasing n. Therefore, if qi lies inside Dq, then qn lies 
inside Dq for all n. Also note that i„ always lies below L, meaning that Dq covers t„. □ 

Lemma 7 Let Dq and Dq be two disks tangent to each other and to an edge segment S at its endpoints. If 
the covering Dq, D\, j = 0, 1,2, . . produced by the edge cover algorithm covers all points of the triangular 
gap T delimited by Dq, Dq and S, then there exists a partition ofT into pieces Ti, « = 0, 1, 2, . . . , such that: 

1. Piece Ti contains Dq: Ti D Dq 

2. Piece Ti is contained inside D\: Ti C D\ 

3. The pieces Ti collectively cover T: UiTi 3 T 
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Proof: Start by assigning points uniquely covered to the only piece that covers it: Ti = TO {D\ — Ui^jD{). 
Next grow each Tj at a uniform rate from their boundaries, but do not permit growth beyond the out- 
circle boundary. Growth of each set is only permitted to consume so-far unassigned points; once a point is 
assigned, it is off-limits for growth. Then Ti,T2, . . ., is a partition of T. □ 



6 Pentagon 

A pentagon has ji — 1/ cos(7r/5) w 1.23607. The lower bound provided by Lemma |21 is -fe ~ 1.11803 (see 
Table ITT^ . Figure [2t^ shows a partition that achieves je, therefore it is optimal. 




We start with the pentagon's inscribed circle Dq and inflate it by to obtain Dl. In each corner of the 
pentagon we nestle five largest possible disks Dq and inflate each by jg to obtain Dl. We choose to make 
Dq and Dl touch each other at the intersection with the corner's bisector, so as to create two symmetrical 
gaps on each side of D^. Cover each of the uncovered edge segments using the edge cover algorithm. The 
algorithm uses 12 iteration steps; therefore, the number of partition pieces is 20476 = 5 * (2^^ — 1) -|- 1, 
the second term counting the big central piece. It is easy to verify that Dq covers the intersection point 
between D^ and Dl; therefore, conform Lemma El the algorithm covers all points interior to the pentagon. 



7 Hexagon and beyond 

A hexagon has 71 = 1/ cos(7r/6) « 1.1547. The lower bound provided by Lemma[21is 79 ~ 1.07735 (see 
Table which is below the critical value 7* of Theorem O Intuitively, this means that it is difficult, 
if not impossible, to achieve 79 for fc-gons for any A; > 6. We use the edge cover algorithm described in 
Section lOl to construct partitions of fc-gons, fc > 6, and compute the best 7 that can be achieved using this 
algorithm. 

For a fixed 7, we partition a fc-gon into pieces with ratio 7 as follows. As before, we start with the 
/c-gon's inscribed disk D^ and inflate it by 7 to obtain D^. In each corner of the fc-gon we place the largest 
possible indisk D^ and inflate it by 7 to obtain We displace Dl along the corner's bisector just enough 
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to capture the corner, as shown in Figure IIUI In this way we create two symmetrical triangular gaps on 
each side of Dq, for a total of 2k triangular gaps that remain to be covered. Cover each such triangular gap 
using the algorithm from section 15. II 




Figure 10: Disk nestled in one corner of the fc-gon. 



We now show how to compute the best balancing 7 for this particular covering. Without loss of generality, 
we consider a k-gon with unit radius indisk and a coordinate system set with the origin at the left corner 
of the bottom horizontal edge. Let 9 = 7r/2 — 7r/fc denote half of the fc-gon's angle. We need to know where 
Di, the inflated central circle, cuts the x-axis closer to origin: 



60 = cot{0) - v/t'-I (16) 
Next, we need to compute the corner indisk Dq-. 

The indisk Dq is tangent to the a;-axis at point 

ai = n cot{0) (18) 
From this, we can compute the point 61 where Dl intersects the x-axis, closer to origin: 

bi — 2ri7 cos(6') 

The edge segment [61,60] is covered using the algorithm from Section Based on Lemma El the gap is 
fully covered if the indisk Dq centered at point (02, r2), with 

02 = {bo + bi)/2 
r2 = (ai - a2)2/4ri, 

covers the apex of the triangular gap. Note that the initial scaled gap value used in equation ((HJ is 
Qq = (60 — ai)/ri. Substituting the expressions for 60 ri (|T7|l and oi ((TH|l . this expands to 

a'o = |;^^!|||^j (cot(g) - Vt^) - com 
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We now solve for 7 that satisfies 



F(ao) < flg 

^'(«o) >0 
(02 - t^f + (r2 - ty) 



(19) 



where {tx,ty) is the intersection point between outcircles D'^ and Dj. Conform Theorem [SJ the first two 
inequahties in H19|l ensure that the edge cover algorithm terminates in a finite number of steps. Conform 
Lemma ini the third inequality in ifT^ ensures that the algorithm covers the entire triangular gap. 




5 10 15 20 25 , 30 

k 

Figure 11: (a) Optimal ratio 7* and single-angle lower bound 70 versus number of vertices k of regular 
/c-gons. 

Solving (|19|l for 7 and k — 6,7 and 8 yields the ratio values shown in Table 11.21 The top curve in 
Figure ITTI shows how the ratio 7* that satisfies (I19II varies with k. The bottom curve represents the single- 
angle lower bound ratio 7^, which is best any algorithm could achieve. As is clear from Figure [TTl the ratio 
achieved by our algorithm is close to the optimal. 

8 Discussion 

We leave open the question of whether optimal paritions can be achieved for fc > 6 with a finite number of 
pieces. 
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